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Abstract 



^ ■ Nonlinear localized excitations in one-dimensional diatomic lattices with cu- 

bic and quartic nonlinearity are considered analytically by a quasi-discreteness 
'^ , approach. The criteria for the occurence of asymmetric gap solitons (with vi- 

brating frequency lying in the gap of phonon bands) and small-amplitude, 

Q ' asymmetric intrinsic localized modes (with the vibrating frequency being 

above all the phonon bands) are obtained explicitly based on the modula- 
tional instabilities of corresponding linear lattice plane waves. The expres- 
sions of particle displacement for all these nonlinear localized excitations are 



O \ also given. The result is applied to standard two-body potentials of the Toda, 

^:^ I Born-Mayer-Coulomb, Lennard-Jones, and Morse type. The comparison with 

previous numerical study of the anharmonic gap modes in diatomic lattices 



for the standard two-body potentials is made and good agreement is found. 
PACS numbers: 46.10+z, 63.20. Pw, 63.20. Ry 



Typeset using REVTeK 



I. INTRODUCTION 

The study of the dynamics of nonhnear lattices and related solitonic excitations has 
greatly influenced by the pioneering works of Fermi, Pasta, and Ulam,lil and of Zabusky 
and Kruskal.B Most of the work in this area has focused on models of one-dimensional (ID) 
monatomic chains with simple interatomic potentials of polynomialETQ, which can approx- 
imate any realistic potential near the equilibrium separation distance of two atoms. This 
description, usually done in a continuum limit, is only valid for a zone-boundary phonon 
mode, i. e. for q, the wavenumber of lattice waves, being near zero or n/do, where d^ is 
lattice spacing. In 1972, TsuruiB proposed an analytical method for studying the nonlin- 
ear excitations of lattices valid in the whole Brillouin zone (BZ). Later on this approach was 
extended by RemoissenelB and Huang.ll2lO Exact analytical solutions for the nonlinear local- 
ized excitations in ID monatomic lattices can be obtained only for the TodaEj and Ablowitz 
and Ladikta lattices, which are discrete completely integrable systems. 

In recent years, the interest in localized excitations in nonlinear lattices has been re- 
newed due to the identification of a new type of anharmonic localized modes.Ej'EJtJ These 
modes, called the intrinsic localized modes (ILM's),c3 or the discrete breathers,B3 are the 
discrete analog of the envelope (or breather) solitons with their spatial extension being 
only of a few lattice spacing and the vibrating frequency lying above the upper cutoff 
of phonon bands.Ej Experimentally, the ILM's have been observed in coupled pendulum 
latticesE^ and electrical lattices.B3 The quantum mechanical aspects of the ILM's have also 
been considered. cItEj However, examination of the ID lattices with standard Toda, Born- 
Mayer-Coulumb, Lennard- Jones, and Morse two-body interatomic potentials demonstrates 
that the ILM's do not appear above the top of the plane wave spectrum. The physical reason 
for this is that the cubic nonlinearity in the Taylor expansion of these realistic potentials is 
too strong. One of the effects of the cubic nonlinearity is that increasing the magnitude of 
the cubic term makes the potentials softer and hence decreases the localized mode frequency. 
The localized mode is destroyed as it approaches the bounding plane wave spectrum.t3'Ey 

Recently, much attention has been paid to the gap solitons in nonlinear diatomic 
lattices.lljcZltlrEJ The concept of the gap solitons was first introduced by Chen and Millsc3 
when investigating the nonlinear optical response of superlattices. For a diatomic lattice, 
the phonon spectrum of the system consists of two branches (acoustic and optical ones), 
induced by mass or force-constant difference of two kinds of particles. Due to nonlinearity 
gap soliton modes may appear as localized excitations with vibrating frequency being in 
the gap of the linear spectrum. Since the gap solitons occur in perfect lattices with discrete 
translational symmetry, a name "anharmonic gap mode" or "intrinsic gap mode (IGM)" was 
given by Sievers and his coUaborators.BtJ'El It is possible that the ILM's and the IGM's may 
be created experimentally in diatomic lattices. Refs. 43 and 44 reported some experimental 
studies of the gap solitons, resonant kinks and the ILM's in a damped and parametrically 
excited ID diatomic pendulum lattices. 

Since the standard two-body potentials of the Toda, Born-Mayer-Coulomb, Lennard- 
Jones, and Morse type have a strong cubic nonlinearity in their Taylor expansion near 
equilibrium position, it is therefore necessary to consider the nonlinear excitations in the di- 



atomic lattices with cubic and quartic anharinoniticity. In their recent contributions, Kiselev 
et alMtJ investigated the anharmonic locahzed modes in ID diatomic lattices with above 
mentioned two-body potentials. By using a rotating wave approximation combined with 
computer simulation, they showed that an ILM does not exist and a nonlinear optical lower 
cutoffgap mode(i. e. IGM) is general feature of these diatomic lattices. Lately, Franchini 
et a/.E3 numerically found that for a potential with the cubic and quartic nonlinearity there 
exists a "critical" K3 (cubic force constant in the potential) value. For small K3, nonlinear 
optical and acoustic upper cutoff localized modes occur, while for large K^ these modes 
disappear and a nonlinear optical lower cutoff mode rises. In a recent work, Bonart et alxB 
investigated the boundary condition effects in the diatomic lattices with cubic and quartic 
anharmonicity. Based on a rotating wave approximation they gave existence criteria for 
the ILM's and IGM's, which are related to the stability properties of linear optical upper 
and lower cutoff phonon modes. These studies posed an interesting problem of how to pro- 
vide an analytical approach which can give not only the explicit criteria for the existence 
of the ILM's and the IGM's as well as other possible nonlinear eacitations for both optical 
and acoustic branches but also the approximate analytical expressions for these nonlinear 
excitations in a unified way. It is this problem that will be addressed here. 

There are several theoretical methods to study the nonlinear localized excitations in 
diatomic lattices (see Ref. 11 and references therein). In this paper we use the quasi- 
discreteness approach (QDA) for diatomic latticeslB to investigate the ILM's and IGM's 
as well as kink-like excitations with small-amplitude in ID diatomic lattices with cubic and 
quartic nonlinear interactions between their nearest- neighbor particles. The paper is orga- 
nized as follows. In Sec. II, the model is introduced and an asymptotic expansion based on 
the QDA is made for the equations of motion. By using the results obtained in Sec. II, in 
Sec. Ill we discuss the solutions of the ILM's and IGM's in a simple and unified way. Some 
explicit criteria and expressions of particle displacement of the ILM's and the IGM's are 
also given in this section. In Sec. IV we apply our results to the standard two-body poten- 
tials from Toda to Morse type and make a comparison with existed numerical experiments. 
Finally, Sec. V contains a discussion and summary of our results. 



II. MODEL AND ASYMPTOTIC EXPANSION 

A. The model 

We consider a ID diatomic lattices with a nearest-neighbor interaction between parti- 
cles. The restriction to the nearest-neighbor interaction is for simplicity and the approach 
can be easily extended to second and higher neighbors. The Hamiltonian of the system is 
given by 



H = E 



1 



-rrii 



dt 



+ V{ui+i - Ui 



where Ui = Ui{t) is the displacement from its equihbrium position of the ith particle with 
mass rrii = m6i^2k + M6i^2k+i {M > m,k is an integer). The potential V{r) is quite general, 
typically it can be the standard two-body potentials of the Toda, Born-Mayer- Coulomb, 
Lennard- Jones and Morse type (for their detailed expressions, see Sec. IV below). We focus 
on displacements with smaller amplitude which can be detected experimentally without in- 
troducing reconstruction or phase transitions in the system. This allows us to Tayler expand 
the potential V{r) at the equilibrium position r = in a power series of the displacements 
to four order .e3 Thus we obtain an approximate K2-K3-K4 potential 

V{r) = ^K^r^ + ^K,r^ + ^-K,r\ (2) 

where K2{> 0), K^ and Ki{> 0) are harmonic, cubic and quartic force constants, respec- 
tively. We assume that the basic features of the weakly nonlinear localized excitations for the 
standard two-body potentials may be obtained by corresponding the K2-K3-K4 potentials. 
Then the Hamiltonian (1) takes the following form 



h=j: 






(3) 



Since each of the standard two-body potentials mentioned above has only one minimum, we 
assume that for the K2-KS-K4 potential (2) there is the constraint 

' <4, (4) 



K2K^ 



unless there are two minima (i. e. double well potential) hence the system may admit some 
types of nonlinear excitations which will not be discussed here. 

If we write U2k = Vn (even particles) and U2k+i = Wn (odd particles), n is the index of 
the nth unit cell with a lattice spacing d = 2do, do is the equilibrium distance between two 
adjacent particles, the system can be split into two sublattices. The equations of motion for 
Vn and Wn are 

d^ 

+K4{Wn - Vnf + {Wn-l - Vnf], (5) 

d^ 
dt^ 

+7^4 [ {Vn - Wnf + {Vn+1 ' Wnf]- (6) 

The linear dispersion relation of Eqs. (5) and (6) is 

Mq) = {h + J2± [{h + J2? - 4/2 J2 sm\qd/2)]-^}-2, (7) 



where I2 = K2/m and J2 = K2/M. The minus (plus) sign corresponds to acoustic (optical) 
mode. At wave number g = the eigen frequency spectrum has a lower cutoff co'_(0) = for 



the acoustic mode and an upper cutoff uj^{0) =00^ = [2(/2 + ^2)]^^^ for the optical mode. 
At g = Tc/d there exists a frequency gap between the upper cutoff of the acoustic branch, 
U-{7i/d) = uji = \/2J2, and the lower cutoff of the optical branch, u+{ir/a) = V2/2. The 



width of the frequency gap is \/2l2 — \^2J2 = V^2!^(l/y^ — l/yM). In linear theory, the 
amplitudes of lattice waves are constants and linear waves can not propagate and will be 
damped when u (the frequency of the waves) lies in the regions uji < uj < uj2 and u > u^. 
Accordingly, these regions are the "forbidden bands" of the linear waves. This property 
of the eigen frequency spectrum results from the discreteness of the system (^. e. discrete 
translational symmery). However, when the nonlinearity in Eqs. (5) and (6) is considered, 
the above conclusions are no longer valid. As localized excitations, some nonlinear modes 
may appear, whose oscillatory frequencies can lie in these forbidden bands of the phonon 
spectrum. 



B. Asymptotic expansion 

We use the QDA for diatomic lattices developed in Ref. 11 to investigate the effects of 
nonlinearity and discreteness of the system. In this treatment one sets 

Unit) = e^S + e\^l + e\^^l + ■■■ (8) 

where e is a smallness and ordering parameter denoting the relative amplitude of the ex- 
citation and ul^l = u^^^^n, T', <Pn) ■ ^n = ^{na — At) and r = e^t are two multiple-scales 
variables (slow variables). A is a parameter to be determined by a solvability condition. The 
"fast" variable, 0„ = qnd — u{q)t, representing the phase of carrier wave, is taken to be 
completely discrete. Substituting (8) into (5) and (6) and comparing the power of e, we 
obtain a hierarchy of equations about v^^ and w^\ {j = 1, 2, 3, ■ ■ ■): 

^' ^a - h (wlil + wl:l_, - 2vl^l) = ^S (9) 



df" 



with 



Mii = 0, (10) 

^ni = 2At^<1 - I^d^^l-i + ^3«l - <i)' - ^sKLi - <!)', (11) 

dtd^n 'din 

^"'" " ^9t9a "'" [ dtdr + ^ da) "'" ' I "^9^ "'""^ + 2! Se "'""V 

+2h{w'^l - <i)(<i - t;2) - 2/3(^gLi - <l)(^S-i - <l - rf J-^S-i) 

d^n 

+/4[«) - vl^lf + (<i_i - v^l)% (12) 



and 






^a - Mvl^l + t^S+i - 2w^£j = Nl^l 



(13) 



with 



N'^l = 0, (14) 

iVS = 2\^y^l + J,d-^vi!l_,, - -^a^i - <n)' + ^3(t^S+i - <],)', (15) 

_2 7or7;(i) -7/;(i)V7;(2) _7„(2)^ _L2 7,r7/^^ _ ,„(i) V./^) _ .,,(2) , ^^„,{i) ^| 

+ J4[«i - <n)' + (^i!Ul - <l)'], (16) 



which can be solved order by order. In Eqs. (9)-(16), we have defined /, = Ki/m and 
Ji = Ki/M{i = 2, 3, 4). The expressions of M^^l and N^^HiJ = 4, 5, ■ ■ ■) need not be 
written down exphcitly here. 



C. Amplitude equations for acoustic and optical modes 

In order to avoid possible divergence for zone-boundary phonon modes, we solve the 
acoustic and optical modes separetely. First we consider the low-frequency acoustic mode 
of the system. For this we rewrite Eqs. (9) and (13) in the form 



Lwl^l = MM^l + Mi/l^^ 






2/2 Ni% 



d^ 



ij) I ,„(i) 



^ + 2h\vl^l = h{wl^l + 



w-A + m£, 



(17) 



where the operator L is defined by 



f,,0) - (— + 2Tn\ (— + 2 Tn\ 7/(^) - h h (v^^^ + 7/^^'^ + 27/(^) 

^'^n,n — I 0,2 + ^-'2 I I 0,2 ^ 2 I "n,ra -'2<>'2 \^"n,n-l + "n,ra+l + ^"n,« 



(19) 



with u^^\{j = 1; 2, 3, ■ ■ ■) a set of arbitrary functions. For j = 1 it is easy to get 



w. 



(1) 



Fio(r,en) + [Fii(r,Ue^'^"+c.c. 



,{i) 



^n:; = ^io(r,^„) - 



/2(l+e- 



-iqd\ 



ujt - 2/2 



-i^ii(r,en)e*^" +C.C. 



(20) 
(21) 



with 0" = qnd — uj^{q)t. uj^{q) has been given in Eq. (7) with a minus sign. The amph- 
tude (or envelope) functions Fiq and Fn are yet to be determined. Fiq is a real function 
representing the "direct current (dc)" part relative to the fast variable 0~ and Fn is a 
complex amplitude of the "alternating current (ac)" part. If K^ = 0, the dc part (Fiq) 
vanishes in this order. For j = 2 (the second order) a solvability condition determines 
A = V~ = dti;_/dg {i. e. the group velocity of the carrier waves) thus ^„ = ^~ = e^nd — Vg't). 
In the third order (j = 3), solvability conditions yields the evolution equations for Fio and 

The detailed expressions of the coefficients in Eqs. (22) and (23) are given in Appendix A. 

Second, we study the high-frequency optical mode excitations. In this case we recast 
Eqs. (9) and (13) into the form 

Lv^l = MN^l + A^Li) + f S + 2^2) Mi% (24) 



^ + 2/2) wiil = J2 {vl^l + v^^^,) + N^^l (25) 



By the same procedure of solving the acoustic mode given above, we obtain 

<l = G'io(r,a) + [G'n(r,a)e*'^" +C.C.], (26) 

-jJl + Qnd\ 



w, 



n,n — Gio{t,^„ 



ui - 2 Jo 



-G'ii(r,a)e*^"+c.c 



(27) 



with 0+ = qnd — u+{q)t. The evolution equations for the amplitudes Gio (dc part of the 
optical mode) and Gn (complex amphtude for ac part of the optical mode) are given by 

^TtGu + 7^c,+ j—j—Gu + P+G^^--G^o + 7+ IGupGn = 0, (28) 

5^^^G,, + a+-^ |Gnp = 0, (29) 

where S,n = ^{nd — K^t) with A = V^ = duj^/dq. The coefficients in Eqs. (28) and (29) are 
also given in the Appendix A. 

Under the transformation 

Fio = (l/e)(?_, Fn = (1/e)/-, (30) 

Gio = (l/e)(?+, Gn = (l/e)/+, (31) 



0, 




(32) 
(33) 




= nd - 


- VgH and 



the nonlinear amplitude equations (22), (23), (28) and (29) can be written in the unified 
form 

r d d , "9 I . ,2 „ 

when returning to the original variables. In Eqs. (32) and (33), 
plus (minus) sign corresponds to the optical (acoustic) mode, respectively. 

Finally, we consider the acoustic mode at g = 0. Noting that Eq. (32) for minus sign is 
invalid at g = for the description of nonlinear excitations since P-\q=o = 1-\q=o = and 
a-\q=o = cxD. This breakdown is due to the fact that at g = an acoustic mode excitation is 
a long wavelength one. In this case a discrete long-wave approximationliJ should be applied. 
By using the same technique used in Ref. 11, for the acoutic mode at g = we obtain a long 
wavelength amplitude equation 

^ + Pu^u + Qu'^u + H-^u = 0, (34) 

where u = dAo/dxn, Xn = nd — ct with c^ = K2d'^/[2{M + m)]. Aq is the leading order 
approximation of f„ and Wn- Since co'-(O) = 0, for the long wavelength acoustic mode 
there is no carrier wave. Thus the displacement of the lattice is purely a "direct current" . 
Eq. (34) without the second term P du/dxn is standard modified Korteweg-de Vries (MKdV) 
equation. Thus (34) is a modified MKdV (MMKdV) equation. Its coefficients are given by 

g^^^f^^y, (36) 



^ d^ ( 2K, 
ri = — 



16 KM + m 



1 mM 



3 {M + mf 



(37) 



III. ASYMMETRIC GAP SOLITONS, KINKS AND INTRINSIC LOCALIZED 

MODES 

When deriving the nonlinear amplitude equations (32)- (34) we have not used any kind 
of decoupling ansatz for the motion of two kinds of particles with different mass. This is one 
of the advantages of the QDA. On the other hand, the nonlinear amplitude equations, which 
are the reduced forms of the original equations of motion (5) and (6) for small-amplitude 
excitations, are valid in the whole BZ (— vr/rf < g < ir/d) except a zero- dispersion point for 
the optical phonon branch.^ Thus one can obtain the gap solitons and ILM's as well as some 
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possible new nonlinear excitations by solving the cutoff modes of the system in a simple and 
unified way. 

(1). Optical upper cutoff mode. For the optical mode at g = 0, we have ujj^ = 
Us = [2K2{l/m + l/M)]i/2, v;+ = 0, x+ = nrf = x„, a+ = -K2(f/[2{M + m)ujsl P+ = 
-K3u;3d/{2K2), 7+ = -3K4^3(l+m/M)V(2i^2), S+ = Kl£/{Mm), and a+ = AK2Ksd{l + 
m/M)'^/{Mm). Eqs. (32) and (33) with plus sign take the form 

4/. + ^«,^/.+7.l/.l7.^0, (38) 

^-9^ = -^\U\' + C„ (39) 



where /+ = /+ exp(— i/S+Cit) with Ci an integration constant and 

7+ 




(40) 



Eq. (38) is standard nonlinear Schrodinger(NLS) equation. It has an uniform vibrating 
solution 

/+ = /oexp(i7+|/o|^t), 

where /o is any complex constant, which corresponds to the linear optical upper cutoff 
phonon mode with a simply frequency shift 7+I/0P and is a fixed point of the system. Note 
that it is possible to eliminate the time dependence by a simply transformation, justfying 
our use of the term "fixed point" for the uniform vibrating solution. In fact, the fixed point 
may also be written as 

/+ = /oexp[2(7+/ot + 0)], 

where /o in this case is any real constant and < < 2tt. In this sense there is a ring 
of fixed point characterized by the different values of the phase, 0. It is easy to show 
that, since a+ < 0, when 7+ < the fixed point is unstable by a long wavelength small 
perturbation. This kind of instability is due to a sideband modulation of the linear optical 
upper cutoff mode. The modulational instability for waves is called the Benjamin-Feir(BF) 
instability^. It is similar to the Eckhaus instability for patterns in extended dissipative 
systems out of equilibrium.^ The BF instability in discrete lattices and related formation 
of solitonlike localized states were already discussed by Kivshar and Peyrard.E3 By this 
mechanism (usually called the Benjamin- Feir resonance mechanismcj) a linear optical upper 
cutoff excitation will bifurcate, grow exponentially at first and then saturate due to the 
nonlinearity of the system. At later stage, a nonlinear localized mode — optical upper 
cutoff soliton is formed. In fact, for 7+ < 0, e. g., 
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-^ < -, (41) 

K2K4 4' ^ ^ 

Eq. (38) admit the envelope (breather) sohton solution 

-r^J r7osech[r7o(a;„-x„„)]exp[-i-|a+|r7ot-i0o], (42) 

where r]Q, 0o and x„(, = norf are constants, uq is an arbitrary integer. Inequality (41) is just 
the condition of the modulational instability for the linear optical upper cutoff mode. From 
Eq. (39) we obtain 

(T_[_Q!_|_ 1 r / \n / \ 

9+ = - . ~ ^0 tanh r7o(a;„ - x„J , (43) 

where the integration constant Ci has chosen as zero as it corresponds to a constant dis- 
placement for all particles. In leading approximation the lattice configuration takes the 
form 

Vnit) = — r^^otanh[r7o(n - no)d\ 

?7osech [r]o{n - no)d] cos{n^st + <Po), (44) 

\- / 

Wn{t) = — T^-z^r]otanh[r]o{n - no)d\ 

0+1+ 

-2— ( -r^ J r^osech [r]Q{n - no)d] cos{Qsst + 0o), (45) 

with 

rt3s = uj3 + -\a+\r]l, (46) 

i. e. the vibrating frequency of the localized mode is above the spectrum of the linear optical 
mode thus above the all phonon bands. Hence (44) and (45) represent an ILM accompanied 
by an asymmetric dc displacement due to the cubic anharmonicity of the system. We call 
it the small-amplitude asymmetric intrinsic localized mode. 

From (44) and (45) we can see that the free parameter rjQ can be taken as an expansion 
parameter, i. e. r]Q = 0(e). By (46) we have r]Q = [2(^^35 — ciJ3)/|a+|]^/^. Thus in our ap- 
proach, the expansion parameter e, used in (8), is propotional to the squareroot of frequency 
difference between the nonlinear localized mode and the linear cutoff phonon mode. 

When 7+ > 0, e.g., the inequality (41) takes opposite sign, the uniform vibrating 
solution of the NLS equation (38) is neutral stable. In this case Eq. (38) admits the dark 
soliton solution 
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/+ = ( -r^ J r7otanh[77o(a;„ - x^o)] exp[z|a+|r7^t - i(j)o]. (47) 

From Eq. (39) we can obtain g^ by integration. In this case we chose Ci in such a waya that 
{dg-^-/dxn)\\xn\=oo = 0. Then we have Ci = o"+|a+|?7Q/((5+7+). Hence we have 

g+ = -^-z — r/otanh ?7o(x„ - x„J . (48) 

The lattice displacement in this case takes the form 

Vn(t) = ^ ^^ ?7otanh[r7o(n - no)d] 

1 

CIq, I \ 2 
-T^ r]otanh[r]o{n-no)d\cos{Q3kt + 4>o), (49) 

7+ / 

u^n{t) = -^-z — r^otanh ?7o(n - no)d\ 

^ ?7o tanh[r7o(n - n^d\ cos(fi3fet + 0o) (50) 



M \ 7+ 



with 



^hk = ^3 r r^Q. (51) 

7+ 

Since 7+ < 0, the vibrating frequency of the kink mode denoted by the expressions (49) and 
(50) is greater than lj^. This is an example of a kink with the vibrating frequency above the 
all phonon bands due to the cubic nonlinearity of the system. 

(2). Optical lower cutoff mode. For the optical mode at g = vr/rf (zone- 
boundary optical phonon mode), one has uj+ = UJ2 = J2K2/m, V^ = 0, a;+ = Xn, 

a+ = K2<f/[2uj2{M - m)], {3+ = -K3u;2d/{2K2), 7+ = -3/^4^^2/(2/^2), h = Kid'' /{Mm) 
and 0"+ = 4:K2K3d/{Mm). In this case we have 

7.-^c..f^-^V (52) 

^^ K2 \K2K4 A ^ ' 



If 7+ > 0, I. e. 



K2K^ ^ 
Eqs. (32) and (33) for plus sign have the solution 
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K'^ 3 

' >i (53) 



(Q, \ 2 ]^ 

-ri J ?7osech [r]o{xn - a;„J] exp[i-a+r]lt - ?0o], (54) 

9+ = - r ~ r]otanh[r]o{Xn - x„J . (55) 

The lattice configuration takes the form 

Vn(t) = — ^^^r]otanh[i]o{n - no)d] 

r^osech [■r]Q{n - no)d] cos{^2st + 0o), (56) 




o" I a I 



Wn{t) = - , ^ r]otanh[r]o{n - no)d] (57) 



with 



1 .2 



^2s = UJ2- -a+%, (58) 

lying in the frequency gap of the phonon spectra between the acoustic and the optical modes. 
It is a typical asymmetric nonlinear gap mode, existing in the diatomic lattices when the 
condition (53) is satisfied. We note that for this mode the displacement of the heavy particles 
only has a kink-like asymmetric dc part. But the displacement of the light particles, besides 
the same type of dc part, has an additional "staggered" vibrational part (z. e. "staggered" 
envelope soliton). We call it the asymmetric optical lower cutoff' gap soliton. The vibrating 
frequency Q2s has the parabola relation with respect to the wave amplitude, denoted by the 
parameter rjQ. The formation of the asymmetric nonlinear gap mode is also the conclusion of 
the BF instability for the corresponding linear optical lower cutoff phonon mode. A further 
discussion for such nonlinear modes in the realistic potentials is given in the next section. 

(3). Acoustic upper cutoff mode. For the acoustic mode at g = tt/c? (zone- 
boundary acoustic phonon mode), one has c^_ = c^i = J2K2/M, V~ = 0, x~ = Xn, «- = 
-K2dy[2uJi{M -m)] < 0, (3_ = -K^uid/{2K2), 7- = -37^4^1/(2X2), 5. = Kld'^/{Mm) 
and a^ = AK2K^d/{Mm). Similarly one can obtain the equations like (38) and (39) with 
5'+, /+, 0-+, (5+, a+ and 7+ changed by g^, f-,a^,6-, «_ and 7_. Here 

a. 2K, ( Kl 3\ 

7- = 7- - P-^- = -r^^i -rr-rr - 7 • (59) 



If 7_ < 0, when 7„ < 0, i. e. 



5. K. '{ KoK, 



Ki 3 



'3 



K.K^ 4 



< 7' (60) 



due to a BF instability of the corresponding linear upper cutoff acoustic mode an asymmetric 
acoustic upper cutoff gap soliton appears, with the vibrating frequency being in the gap of 
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the phonon spectra between the acoustic and optical modes. Otherwise an acoustic upper 
cutoff kink vibrational mode occurs. We can readily written down the explicit expression of 
the lattice configuration in this case, but for saving space it is omitted here. 

(4). Acoustic lower cutoff mode. This is a long wavelength mode without any 
carrier wave, because when g = we have uj^ = thus (f)~ = 0. The lattice displacement 
only has dc part and its evolution is controled by the MMKdV equation, given by (34). A 
single-soliton solution of Eq. (34) is given by 

y^ = ^ (Ql) 

P + y/P'^ + 2Ak'QH cosh [2K{nd - (c + 2AK'^H)t)\ ' ^ ' 

where k is an arbitrary constant. Making the transformation u = —P/{2Q) + U{(^n,t) with 
Cn = Xn + PH/{AQ), Eq. (34) becomes 

It is the standard MKdV equation and can be solved by the inverse scattering transform. 
For the explict expressions of the kink and breather solutions of the MKdV equation we 
refer to Ref. 11. 

Needless to say, in addition to the cutoff modes considered above, our approach devel- 
oped in Sec. II can also be used to discuss the nonlinear localized excitations for g 7^ and 
q 7^ vr/d {%. e. the intra-band modes), which will be done elsewhere. 



IV. APPLICATION TO REALISTIC TWO-BODY NEAREST-NEIGHBOR 

POTENTIALS 

In this section, we apply the general results obtained above to the diatomic lattices 
with the standard two-body nearest-neighbor potentials to see whether the anharmonic gap 
modes and the ILM's can appear or not. This can be easily done by using the existence 
criteria given by (41) (for ILM's), (|53| ) (for optical lower cutoff gap solitons), and (pO| ) (for 
acoustic upper cutoff gap solitons). Four standard interatomic potentials area 



1. Toda: 



V(r) = -e-""- + ar--, (63) 

b b 



where a and b are coefficients, such that a6 > 0. r is the deviation of the relative 
interparticle separation from its zero-temperature equilibrium position. 



2. Born- Mayer- Coulomb: 

"0 



d\ 



2 



r -\- do 
13 



V+Poe-^/''°+rfo- 



(64) 



where um is the Madelung constant, q is the effective charge, do is the zero-temperature 
equihbrium distance between adjacent particles, and po is the constant describing the 
repulsion between atoms. 



3. Lennard- J ones: 

V{r) = a 



, \ 12 / , \6 

"0 \ o / "0 \ 



r + doj \r + do 



+ 1 



(65) 



where a is the constant determing the potential strength. 
4. Morse: 

\/(r) = p(e-""-l)', (66) 

where P and a are constants determing the strength and the curvature of the potential, 
respectively. 

For weakly nonlinear excitations we can Taylor expand these potentials at their equi- 
librium position (r = 0) to obtain the force constants. They are defined by 

1 /d^V\ 

with j = 2, 3, 4, ■ ■ -. Thus we have the values of K2, K^, K4 and i^|/(i^2-^4) which are given 
in Table 1. Notice that K^ < for these potentials except that the Toda one with a and b 
being both negative (note that we require that K2 and K4 are positive). The parameter / 
for the Born-Mayer-Coulomb potential in Table 1 is defined by 



(dl - Qpl) 



{do-2po){dl-24p: 



(6J 



Generally, we have / > 1. For example"'"'", for KI, (io=3.14A, po = 0.26A, one has 
/ = 1.1171; for KBr, do = 3.29A, po = 0.334A, we have / = 1.1328; for Lil, do = 3.0A, 
Po = O.374A, one has / = 1.1487. By these results as well as the criteria given by (41), (53) 
and (60), we make the following conclusions: 

(1). The condition (53) is satisfied by all these standard two-body potentials. Thus 
the asymmetric optical lower cutoff gap soliton mode given in the expressions (56) and (57) 
does exist in the diatomic lattices with the Toda, Born-Mayer-Coulomb, Lennard- Jones and 
Morse type interatomic interactions. This result agrees with the conclusion by numerical 
study presented by Kiselev et alMM 

Shown in Fig. 1 is the dc (denoted by the solid circles) and ac (denoted by the open 
circles) amplitude patterns of light particles for the optical lower cutoff gap soliton mode. 
The amplitude pattern for the heavy particles, which is not shown here, only has dc part 
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similar to that of the dc part for the hght particles. The Born-Mayer-Coulomb potential 
for KBr-like parameters is used with m/M = 39/80. Comparing the panel (a) of Fig. 1 in 
Ref. 34 with our result shown in Fig. 1 here, we see that the lattice configuration obtained 
analytically by our QDA is basically the same as the corresponding result given by the 
rotating wave approximation plus computer simulation, used by Kiselev et a/.o Furthermore, 
the parabola relation between the vibrating frequency and the amplitude of the nonlinear 
optical lower cutoff gap mode, given by (58), is in accordance well with the numerical 
results (see the Fig. 3 of Ref. 34). 

Shown in Fig. 2 is the amplitude (maximum absolute value) ratio (|y4rfc|/|^ac|) of the dc 
part (I Add) to ac part (l^ad) of the light particle displacement for the optical lower cutoff 
gap soliton mode. We see that the amplitude ratio is the function of the mass ratio m/M. 
The larger is m/M, the larger is the amplitude ratio. The amplitude ratio grows very fast 
as the mass ratio approaches 1. 

(2). Notice that the existence criterion of an ILM (the optical upper cutoff soliton) is 
the inequality (41), i. e. the necessary condition for the occurence of the ILM is 

K'^ 3 

^' <-. (69) 



K2K^ 4 

From Table 1 we see that in general i^|/ (7^2-^4) > 1- Thus for all the standard two-body 
potentials from Toda to Morse type, an ILM is impossible. This conclusion concides also 
with the numerical results of Refs. 33 and 34. 

(3). An acoustic upper cutoff gap sohton is not possible for all these standard two- 
body potentials, because its existence condition (60) (the same as (41) ) can not be satisfied. 
However, acoutic and optical upper cutoff vibrating kinks are possible excitations for these 
diatomic lattice systems. 

(4). Nonlinear acoustic lower cutoff excitations are long wavelength modes and are 
governed by the MMKdV equation (34). Since Q and H are positive and P is negative (due 
to K3 < 0) for all the two-body potentials (except the Toda one with a and b being both 
negative), the soliton amplitude (see (61) ) in the presence of cubic nonlinearity {K3 7^ 0) is 
smaller than the soliton amplitude in the absence of the cubic nonlinearity (i^s = 0). 

(5). For general K2-K^-Ki^ potentials, by the criteria given by (41), (53) and (60), we 
conclude that in weakly nonlinear approximation the "critical" value of K^ / {K2K4/) for the 
transition from an optical upper cutoff kink to an optical upper cutoff soliton (ILM) and for 
the occurence of an optical (acoustic) lower (upper) cutoff gap mode is 3/4, independent of 
m/M. Because one of the conditions (53) and (60) must be satisfied, we have the conclusion 
that for any nonlinear diatomic lattice, gap solitons always occur. Our analytical results 
support the numerical findings of Franchini et a/.H for small cubic nonlinearity, nonlinear 
optical and acoustic upper cutoff localized modes appear, while for large cubic nonlinearity 
a nonlinear optical lower cutoff mode rises. When K-^ = O^the corresponding theoretical 
and numerical results for harmonic plus quartic potentialstiJo are recovered. 
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It is interesting to note that the criterion (41) and (60) for the occurence of the asymmet- 
ric optical and acoustic upper cutoff sohtons in the diatomic lattices are the same as that for 
the appearence of the upper cutoff solitons in monatomic lattices, given by TsuruiH (Eq. (4.8), 
ujI = 4), Flytzanis et alE (Eq. (5.5), kD = tt) and Flachll(Eq. (3.24), v^ = 0). This criterion 
was also given implicitly in Ref. 10 since the envelope soliton solution (23) in Ref. 10 requires 
sgn(PQ) > 0. 

In addition, our results show analytically that for the potentials with cubic and quartic 
nonlinearity an nonlinear localized excitation always consists of dc and ac parts. The ap- 
pearrence of the dc part is a direct conclusion of the asymmetry in the potentials. This fact 
is already known by using different approaches (see e. g. Refs. 33 and 34). 

Recently, Bonart, Rossler and Pagec3 considered the boundary condition effects in the 
diatomic lattice with cubic and quartic anharmonicity and gave existence criteria for the 
ILM's and the IGM's not restricted to small amplitudes by using a different approach. By 
comparison we can see that our instability threshold criterion for the linear optical upper 
cutoff phonon modes is in accordance with their corresponding one [Eq. (10) in Ref. 45)] 
when exact to A^^ [remenber that in our notation Aac is the amplitude of the ac part of the 
nonlinear excitation, e.g. in (44) we have Aac = '2{ci^+/7+Y^'^Vo\- It is easy to show that 
our frequency formulas for the nonlinear optical (including upper and lower) cutoff modes, 
i. e. Eqs. (46) and (58), also coincide with the corresponding ones in Ref. 45 [Eqs. (9) and 
(12)] if the terms propotional to A^^ are neglected. Our results support also the criterion 
(14) in Ref . 45 for the appearance of nonlinear localized modes. In addition, our analytical 
approach based on QDA may provide much more insights such as for the occurance of 
nonlinear localized excitations in acoustic branch although it has the limitation of small 
amplitude. 



V. DISCUSSION AND SUMMARY 

Based on the QDA for diatomic lattices we have analytically studied the nonlinear lo- 
calized excitations with small amplitude in the diatomic lattices with cubic and quartic 
nonlinearity. The results are quite general and allows us directly to obtain many different 
types of nonlinear excitations in a unified way. Starting from the nonlinear amplitude equa- 
tions given in (32) and (33), the existence criteria for the optical upper cutoff solitons and 
asymmetric gap soliton modes have been explicitly provided in (41), (53) and (60). The an- 
alytical expressions of particle displacements for all these nonlinear localized modes are also 
given. The theoretical results have been applied to the standard two-body nearest-neighbor 
potentials from Toda to Morse type and agreements well with the previous numerical findings 
have also been found. 

Most of the existed analytical studies for the nonlinear excitations in diatomic lat- 
tices involved so-called "decoupling ansatz" , in which some relations were assumed between 
the displacement of light particles and that of heavy ones before solving the equations of 
motion (see Ref. 11 and references therein). Pnevmatikos et a/.H investigated the soliton 
dynamics of nonlinear diatomic lattices by using the decoupling ansatz. For long wave- 
length {i. e. q = 0) excitations this can be done without big difficulty. But for envelope-type 
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excitations and in the case of cubic nonlinearity, the concrete form of the decouphng ansatz 
is not easy to determine and the analytical calculation for the coefficients in amplitude equa- 
tions, like (32) and (33), is also heavy. Except several numerical studies, such an analytical 
calculation has never been accomplished. Recently, it was shown that the decoupling ansatz 
is completely unnecessary and can be derived by the QDAllll. In addition, the QDA has 
many other advantages. For example, the results obtained by the QDA, though restricted 
to small amplitudes, are valid in the whole BZ except at the zero-dispersion point of the op- 
tical phonon branch ( see Ref. 46). Thus one can obtain all solutions for nonlinear cutoff and 
non-cutoff modes in a simple and unified way; the method is quite general and can be ap- 
plied to the other lattice systems. An extension based on the QDA for magiietic gap soliton 
excitations in alternating Heisenberg ferromagnets has been given recently.^ In the present 
work, we use the QDA to consider the nonlinear localized excitations in the diatomic lat- 
tices with cubic and quartic nonlinearity. The detailed expressions of the coefficients in the 
amplitude equations (32) and (33) as well as (34) thus various types of nonlinear excitations 
are explicitly given. 

The study of stability by using nonlinear amplitude equations is widely employed in 
pattern formation in systems out of equilibrium.E3'E3 In present work, by a similar idea we 
have studied the stability of the linear optical and acoustic (upper and lower) cutoff phonon 
modes. For the i^2--^3--^4 potential, we have obtained the existence criteria for the occurance 
of linear cutoff phonon mode-related nonlinear localized excitations: If K^ / {K2K4) < 3/4, 
we have the acoustic and optical upper cutoff solitons. Otherwise one has only the optical 
lower cutoff (gap) solitons. The later case happens for the standard two-body potentials 
from Toda to Morse type. The formation of the asymmetric gap solitons in the standard 
two-body potentials is the result of BF modulational instability of the optical lower cutoff 
phonon modes. Our results show that for any nonlinear diatomic lattice a gap soliton always 
occurs. The reason for this is that the curvatures of the acoustic and optical branches of 
the phonon spectra in the vicinity of the BZ edge have different signs. This means that if 
the optical mode at the BZ edge supports a lower cutoff gap soliton, the acoustic mode at 
the BZ edge supports an upper cutoff kink subject to the same sign of the nonlinearity, and 
vice vesa. 

For large-amplitude excitations we should extend the present QDA to a higher-order 
approximation. It is expected that higher-order correction terms, likely |/±|^/± and 
{d/dx^)\f±\^ {as well as some higher-order dispersion terms), will be added respectively 
into Eqs. (32) and (33). Then the existence criteria for the nonlinear localized excitations, 
i. e. (41), (53) and (60), will be modified to be amplitude-dependent. 

The generation of an asymmetric IGM in nonlinear 3D diatomic lattices has been con- 
sidered recentlytZlo. These studies reveal the possibility of observing experimentally the 
IGM's in real (3D) crystals. An macro-analogy of the diatomic lattices with cubic and quar- 
tic nonlinearity is to use a chain of magnetic pendulums. Russell et alM reported moving 
breathers in such a system quite recently, but they did not pay attention to nonlinear cutoff 
modes. If the mass of the magnetic pendulums is arranged in a alternating way, one can 
obtain a diatomic lattice with an asymmetric intersite (dipole-dipole interaction) potential. 
For small amplitude excitations the intersite potential reduces to a K2-KS-K4 one. Thus it 
is possible to observe the asymmetric IGM's by using this macro-diatomic lattice system. 
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APPENDIX A 



The detailed expressions of the coefficients of Eqs. (22) and (23) are given by 



hJ2d^ - 2(J2 + J2) 



l2J2d sin{qd) 

2UJ4I2 + J2 - iO^-) 



a_ = 4/3/3^ 



ul - 2 Jo 



+ 



472/3^ 



ui - 2/0 ui - 2L 



(c(j^ + 2/2 COS qd) 



32/272-^3^? ^- singe? 



a. 



4(12 + J2- Scu^-)iVg-y - d^ljcul - 2l2){ujl - 2J2) - 2/2J2] 



^ _ljuj_d 
^-~ 2/2 ' 

3J4C<J_ 



7- 



2J2 



^UJ-iJt 



cj2 - 2 J2 / 



h-jl 



ui - 2/2 



oo^_ sin^{qd) 



[UJ- 



2I2W2 + J2-UJI) 



-hL^ 



Js{Is + /2L_ COS qd) 

2cu2 _ /„ 



L 



l2.h[cosqd-{2ujl-l2)/J2] 



{2ujI - J2){2ujI -h)-l2J2 cos2 qd' 
The coefficients of Eqs. (28) and (29) are given by 

l2J2dsm{qd) 



S+ = l2J2d' - 2(/2 + J2 

,ujI-21 



2UJ+{I2 + J2 - ^H 



a+ = 4J2 J3(i 



2 , 4/2J3ci ,2,07 j\ , ?>2l2J2hV+uj+smqd 
+ ^^ ^^(w+ + 2 J2 cos qd) + .2 oT^2 ' 



7+ 



Co'.l. — 2^2 '^+ — 2J2 V^+ — ^"^2 

4(/2 + J2 - 34)(v;+)2 _ ^2 [(^2 _ 2j2)(^2 _ 2J2) _ 2/2 J2] 



4a;^('ci;? 



-(2 — ■J'? 



J^u+d 
2J2 ' 

3/4CJ_| 



Cj2 _ 2/2 / 
1 + ^ 1 



c<jjj_ sin^ (qd) 



2/2 



J2 

cji - 2 J2 



[UJi 



2J2){l2 + -h - UJ'\ 



-J,L, 



h{J3 + J2L+ cos qd) 
2ujI - J2 



L 



J2h[cosqd - {2uj\ - J2)/l2 



{2uj\ — l2){2uj\ — J2) — I2J2 cos^ qd'' 
where uj±{q) has been given in Eq. (7). 



(A.l) 

(A.2) 

(A.3) 
(A.4) 



(A.5) 
(A.6) 



(A.7) 
(A.8) 
(A.9) 
(A.IO) 



(A.ll) 
(A.12) 
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TABLE I. The force constants K21 K^, K4 and the value of K'^/[K2K/i) for the stan- 
dard two-body potentials from the Toda, Born-Mayer-Coulomb (B-M-C), Lennard- Jones (L-J), and 
Morse type. 
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FIG. 1. The dc and ac amplitude patterns of light particles in a diatomic lattice for the optical 
lower cutoff gap soliton mode. The solid (open) circles denote dc (ac) part of the lattice displace- 
ment. The Born-Mayer-Coulomb potential for KBr-like parameters is used with m/M = 39/80. 
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FIG. 2. The amplitude(maximum absolute value) ratio (|ylrfc|/|^ac|) of the dc part (|^(;;c|) to 
ac part (|^ac|) of the light particle displacement for the optical lower cutoff gap soliton mode as a 
function of the mass ratio m/M. 
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